Friedel Oscillations and Charge-Density Waves Pinning in Quasi-one-dimensional 

Conductors: An X-ray Access. 
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We present an x-ray diffraction study of the vanadium-doped blue bronze Ko.3(Moo.972 Vo.o2s)03. 
At low temperature, we have observed both an intensity asymmetry of the ±2/cf satellite reflections 
relative to the pure compound, and a profile asymmetry of each satellite reflections. We show that 
the profile asymmetry is due to Friedel oscillation around the V substituant and that the intensity 
asymmetry is related to the charge density wave (CDW) pinning. These two effects, intensity and 
profile asymmetries, give for the first time access to the local properties of CDW in disordered 
systems, including the pinning and even the phase shift of Friedel oscillations. 

PACS numbers: 71.45.Lr, 61.10.-i, 61.72.Dd 



The effects of impurities on the physical properties of 
crystals is one of the major issues of solid state physics. 
In metals, impurities are known to originate Friedel os- 
cillations O (FOs) of the electronic density and to af- 
fect generic ones like charge-density waves (CDWs) 
In quasi-one dimensional (ID) metals both features oc- 
cur. Due to its divergent susceptibility at twice the Fermi 
wave vector 2kp, a ID electron gas is unstable with re- 
spect to the formation of a 2fcp-CDW, which, coupled to 
the lattice, induces a 2/cF-periodic lattice distortion. The 
2fcp-satellite reflections arising from this distortion pro- 
vides a rich variety of information on the CDW. While 
the peak intensity and its position give the amplitude and 
the wave length of the CDW, the satellite reflections pro- 
files is related to the phase-phase correlation function of 
the CDW. In CDW materials like blue bronze K0.3M0O3 
H or transition metal trichalcogenides like NbSe3 [0 , || , 
diffraction experiments have clearly demonstrated the 
loss of the CDW long range order upon doping. But de- 
spite this recent progress in the description of the CDW 
structure, the microscopic nature of the pinning is still 
largely unknown. Recently, two new methods have been 
found to access the local properties of CDW in disor- 
dered systems Q . The first one is the analysis of the 
intensity asymmetry (IA) of the satellite reflections at 
+2kp/— 2kp from a Bragg reflection, which gives infor- 
mation of the value of the phase at an impurity position. 
Another method, that we describe here for the first time, 
consists in analysing the profile asymmetry (PA) of the 
individual peaks. In this letter, we shall demonstrate that 
both effects, intensity and profile asymmetries, can be 
observed simultaneously and that they provide a unique 
tool to access the FO's and the pinning. 

Let us first consider the general case of a modulated 
crystal containing impurities in concentration c. The 
substitution disorder is characterized by cr(r), which is 1 
if an impurity sits in position r and otherwise ^| . The 
displacements of the atoms from their regular positions 
are u(r). The total diffracted intensity at the scattering 



vector Q = G hkl +q close to one of the reciprocal wave 
vector G hk i reads, @: 



J(Q) = I L + I d + I a 

= A/ 2 (a q a_ q )+7 2 (|Q.u 



(1) 

2/A/Im(a_ q Q.u q ) 



where (..) denotes an average on both disorder and ther- 
mal fluctuations. Af = fi — f is the difference be- 
tween the impurity and the host atom form factors, 
/ = cfi + (1 — c)/ is the average form factor and u q 
and <7 q are the Fourier transforms of u(r) and <r(r) — c 
respectively, q, which is any wave vector of the first Bril- 
louin zone, is close to 2kp in the present case. 

The first term gives the scattering intensity due to dis- 
order, usually called the "Laue scattering". The second 
term represents the intensity due to atomic displacements 
i.e. the usual Fourier transform of the correlation func- 
tion (u(O)u(r)). The third term arises from the coupling 
between the disorder and the displacements. Remark- 
ably, I a is linear in u q and odd in q, which changes 
the relative intensities at G/jjy+q and Ghki— q with re- 
spect to that of the pure crystal. An intensity asymme- 
try occurs around the Ghki reciprocal position. Phys- 
ically, this term is due to an interference between the 
reference wave scattered by the impurity and the wave 
scattered by the atomic displacements originated by the 
same impurity. In that sense, this phenomenon can be 
seen as an holography of the lattice distortion around the 
impurity. I a cancels out when there is no coupling be- 
tween the impurity positions and the atomic displace- 
ments ((<T_ q Q.u ) = (cr-q) (Q.u ) = 0). However 
such a coupling does exist for FO's and CDW when it 
is pinned. In a simple model of strong pinning, where 
the lattice distortion has an amplitude U2k F and a phase 
ipo at the impurity site, it has been shown that 

UQ = G hkl ±2k F ) ~ T/A/(Q-u 2kF ) cos </?(,. (2) 

This gives an intensity asymmetry of the +2kp / ~2\ip 
satellite reflections and provides a unique way of cvalu- 
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ating cosc/?o- In organic materials, "white lines" due to 
strong intensity asymmetry have been observed and used 
to determine the phase <po , |J . With the same formal- 
ism, it can be shown that if the phase of the CDW un- 
dergoes a deformation around the defect, a profile asym- 
metry of each satellite is produced . The same effects 
is produced by the notorious phse shifts of FO's. In this 
communication we report on how these ideas have been 
tested in doped blue bronzes, which can be considered, 
with NbSe3 as a prototype of charge-density wave mate- 
rials. 
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FIG. 1. High resolution profiles of the large-fc 

(1, 9 + qb, 6.5) and (1, 11 — §(,, 6.5) peaks at 60 K. The position 
of the q c - and q s -scattering are indicated on the x-axis. The 
vertical arrows indicate the intensity asymmetry of the q c - 
and q s -scattering. 

Quasi-one-dimensional blue bronze K0.3M0O3 under- 
goes a Peierls transition at T p =183K (ljj] . Below T p , 
satellite peaks appear at the reduced wavevector q c , 
with q c = (1,0.748,0.5) at 15K |||. In order to maxi- 
mize the effect of impurities on the scattering, we have 
investigated heavily vanadium-doped blue bronzes (2.8 
at.% V). X-ray diffraction studies were performed with 
a three-circle diffractometer on a 12 kW rotating an- 
ode {XiAoKa =0.711 A) and with a four-circle diffractome- 
ter at the BM2 beamline in ESRF for high resolution 
measurements (A=0.7A and A=0.8A). Previous studies 
had shown that at low temperature, the 2/cF-satellite 
reflections observed at the reduced position q c =(l, 1 — 
qb = 0.68,0.5) were broadened, indicating the lost of 
the CDW long range order j|. At last, the measured 
2k p = (1 — qb)b* (b* = 2Tr/b) was in accordance with the 
2kp value expected from the change in band filling, due 
to the substitution of 2.8% V 5+ for Mo 6 + g. 

In order to gain evidence for an intensity asymmetry, 
we have measured the profiles of 2kp/ — 2kp scattering at 
different reciprocal positions and carefully compared the 
results obtained in doped and pure crystals. Typical high 
resolution scans in the chain (b) direction around two 
satellite reflections are shown in Fig. 1, after background 
substraction. A casual inspection of Fig. 1 clearly reveals 



a profile asymmetry of the satellite reflections, already 
noted in Ref. [§. This asymmetric peaks can be analyzed 
by a sum of two components, a broad one, located at the 
reduced position q s = (1,1 — g s ,0.5), where g s =0.5 at 
15 K, and a narrow one, corresponding to the regular 
q c -scattering by the CDW. 



w 15 - 




-O.B 



-0.6 -0.4 
k ( b* ) 



FIG. 2. Profile of the (13, —1 + qt, —6.5) diffuse scattering 
at 300 K and 15 K. The downwards (upwards) arrows point 
towards the q c -scattering (q s ), respectively. The solid and 
dotted lines are the results of fit described in the text. 

This analysis is confirmed by the temperature varia- 
tion of the satellite profile. Fig. 2. shows that the q s 
broad component is still present at 300 K, with the same 
intensity. The curves in fig. 2 are the result of a fit by a 
sum of Lorentzian-squared functions, convoluted to the 
resolution function. From this analysis, the temperature 
dependences of 1 — qb and 1 — q s have been extracted and 
are displayed in Fig. 3. 




FIG. 3. Temperature dependence of 1 — qt (right scale) and 
1 — q s (left scale) . 

Another striking feature of Fig. 1 is the difference be- 
tween the +2kp and — 2kp profiles. Consistenly with the 
previous analysis, this has been interpreted as due to an 
intensity asymmetry I' a of the q s -scattering, as indicated 
in Fig. 1. At last, a comparison with the satellite reflec- 
tion intensities of the pure blue bronze (not shown here) 
reveals an intensity asymmetry of the q c -scattering [ [l3"| . 
This is evidenced by an intensity ratio I-2k F / hk F 20% 
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lower than that of the pure compound at low temperature 
|^3| . This effet is indicated by I a in fig. 1. 

As far as the low temperature CDW correlation lengths 
are concerned, they can be estimated from the inverse of 
the Half Widths at Half Maximum (HWHM) of the q c - 
scattering. Along the chains, one finds 4. =28.7 A at 15 
K [l^ ]. The ordered domains are quasi-one-dimcnsional 
and they contain on average one V substituant @. 

Let us now present a simple model, which indicates 
that the q s -scattering is due to Friedel oscillations. At 
T = K, the oscillating part of the electronic density of 
a metal at large distance r from an impurity of charge -Z 
located at r = reads [fill : 



5p{r) 



cos(2kpr + 77) 



(3) 



where D is the space dimension and 77 is the phase shift 
of the electronic wave functions at the Fermi level. With 
these definitions, 77 is related to Z through the Friedel 
sum rule Z — -77. 

Physically, this relation means that the additional 
charge -Z has to be screened by the conduction electrons, 
by bringing an opposite charge Z in the vicinity of the 
impurity. In ID, the lattice distortion associated to the 
FO reads : 



u(x) = u a 



exp(- |x| /£) 



sm(2kpx + x( x ))- 



(4) 



where r = |x|. In this expression we have introduced 
a damping length £, due to the temperature and/or the 
disorder, and the phase function x{ x )> whose limits are 
X(±oo) = ±rj, consistently with equation (||). x( x ) 
jumps between these values on a distance roughly equal 
to the extension of the potential of the impurity, i.e. less 
than the lattice spacing. At last, let us remark that a de- 
crease of the distance between neighbors corresponds to 
a decrease (increase) of the hole (electron) density. The 
two modulations are thus in quadrature, in agreement 
with expressions and (Q). 

For the sake of simplicity, we shall model the x-ray 
scattering by taking the ID lattice Fourier transform u q 
of u(x), which can be calculated exactly fl(i|| . This ID 
approximation is justified by the small tranverse correla- 
tion lengths of the q s -scattering. The intensity diffracted 
by the FO, given by the Id ~ |u g u_ g | term, mainly de- 
pends on £ and 77. Solid lines in fig. 4 represent the 
calculated intensity, for 77 = ±tt/2 and the blue bronze 
low temperature value of 2kp = 0.75b*. Indeed, this 
value of r\ corresponds to the screening of a single charge 
impurity, which is the case here as discussed below. For 
large values of £, the intensity exhibits a sharp disconti- 
nuity at ±2fcp and a long tail in the small-q direction. 
These discontinuities are smoothed for smaller values of 
£, which slightly shifts the maxima from the ±2kp po- 
sitions. As expected from the general analysis 0, the 
phase distortion due to the FO phase shift gives rise to 
a profile asymmetry of the satellite reflections. 



In fact, the total intensity diffracted by the FOs in the 
real crystal is given by ([[]). As the Laue scattering is 
weak here, the only additional term is the third one : 
I' a . The dashed line of the lower part of Fig. 4 gives the 
total scattering intensity in the A/ sin 77 < case. As 
expected, an IA of the +q/ — q is obtained, which gives 
the possibility of measuring sin 77. Indeed, using (fij), the 
IA term can be written: 



I' a (Q = G hkl ±q) ~ ±/A/(Q.u )sin77. 



(5) 
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FIG. 4. Intensity of the diffraction by a Friedel oscillation 
for £=1000 b (top) and £= b (bottom). The dashed curve is 
obtained by adding the linear term I' a for A/ sin r\ < 0. 

A comparison of Figs. 1 and 4 strongly suggests the 
q s -scattering to be due to FOs around the single charged 
V impurities. From the position of the q s -scattering with 
respect to 2hp, the damping length can be estimated as 
£ ~ 8 A. The FOs are present at ambient temperature, 
as shown in Fig. 2, and still exist in the CDW phase with 
the same characteristics. Concerning the q c -scattering, 
it corresponds to the regular scattering by the 2k,F lattice 
distortion at larger distance from the impurity. The cor- 
related variation of the q c - and the q s -scattering maxima 
(Fig. 3), both related to the thermal variation of 2kp ||, 
supports this interpretation. 

The observation of intensity asymmetries on both q c - 
and q s -scattering gives evidence of a coherence between 
the position of an impurity, and respectively the FO and 
the CDW. Using eq. (|) and (|) and the signs of I a and I' a 
obtained from the experiment, one can readily calculate 
the signs of costpo and sin 77. Considering that A/ = 
fv ~ I Mo < and Q.u 2kf , > for positive k- values, one 
gets cost^o > and sin 77 > 0, which is consistent ipo ~ 
and 77 ~ 5. Both values correspond to a local increase 
of the hole density around the impurity, which is natural 
for the V substituant. Indeed, the V 5+ atom provide a 
negative charge (Z=l) with respect to the molybdenum 
Mo 6+ background. This charge has to be screened by a 
positive charge, which induces a local increase of the hole 
density as found experimentally. The r\ — \ value is that 
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expected from the Friedel sum rule with Z=l 

Fig. 5 displays a tentative representation of the 
CDW/FO structure around a V atom according to the 
present experimental results. FO's dominate the charge 
modulation in a tiny 2£ ~ 16 A region around the impu- 
rity. Consistently with this result, this damping length 
can be estimated at £ ~ 9 A from ~ t^ 1 +£ ~ 1 , where 
£o = fivp/A ~ 13 A is the electronic coherence length 
fy}, and I = 27 A is the average distance beetween V 
substituants along the chains. At larger distances, the 
CDW dominates the charge oscillation. However, the 
CDW-phase extrapolated at the impurity site is (fo ~ 0, 
as indicated by the dotted line. This description is consis- 
tent with the strong pinning picture as proposed in refs. 
||l6| . In addition, the region of mismatch between the 
CDW and the FO corresponds to the q ~ £ _1 wave vec- 
tors, which indicates that the whole ±2kp profiles con- 
tain also information on this important crossover region. 
Finally, let us point out that these ±2kp- asymmetric 
profiles, which clearly depends on FOs, CDWs and their 
interferences, are of great importance for the understand- 
ing of the sliding of CDWs and open a unique access to 
the properties of CDW at both microscopic (FO) and 
mesoscopic pinning) scales. 



asymmetry analyses of the x-ray scattering allows one to 
study the subtle interplay between structural and elec- 
tronic properties in disordered systems. 

We are indebted to J.-F Berar and R. Moret for their 
help during the synchrotron radiation experiments. 




FIG. 5. Bottom : sketch of the oscillating part of the hole 
charge density around an impurity placed at the origin, as a 
function of the cell index n. Circles represent the Mo atoms. 
Top : Phase \(x — nb) of this charge density, in radian. 

The conclusion of this study is the evidence of Friedel 
oscillations in the vicinity of vanadium atoms in V-doped 
blue bronze. The observation of intensity asymmetries 
provides evidence of a coherence between the impurity 
position and the CDW/FO. Moreover, the observation 
of the profile asymmetry clearly indicates the presence 
of phase distortions, corresponding to a local decrease of 
the electronic density. Additional experiments on crys- 
tal with different doping levels are planned in order to 
elucidate more quantitatively the microscopic features 
of the pinning of CDWs in low-dimensional materials. 
More generally, we have demonstrated how methods of 



[1 

[2. 

[3; 

[4; 

[5; 
[<: 

[8 

[9 
[10 

[11 

[12 
[13 
[14 

[15, 
[16 



J. Friedel, Phil. Mag. 43, 153 (1952); Adv. Phys. 3, 446 
(1954); Nuovo Cimento Suppl. 87, 287 (1958). 
Low-Dimensional Properties of Quasi- One- Dimensional 
Materials, edited by P. Monceau (Reidel, Dordrecht, 
1985). 

S. Girault, A. H. Moudden, J.-P. Pouget, and J.M. Go- 
dard, Phys. Rev. B38, 7980 (1988). 

E. Sweetland, A. C. Finnefrock, W. J. Podulka, M. Sut- 
ton, J. Brock, D. Dicarlo and R. E. Thorne, Phys. Rev. 
Lett. 65, 3165 (1990). D. Dicarlo, R. E. Thorne, E. 
Sweetland, M. Sutton, J. Brock, Phys. Rev. B50, 8288 
(1994) . 

S. Rouziere, S. Ravy, J.-P. Pouget and R. Thorne, Phys. 
Rev. B59, 15121 (1999). 

S. Ravy, J.-P. Pouget and R. Comes, J. Phys. I (France) 
2, 1173 (1992). 

S. Ravy, J.-P. Pouget, J. Phys. IV (France) 3, 109 (1993). 
M. A. Krivoglaz, "Theory of X-ray and Thermal- Neutron 
Scattering by Real Crystals", Plenum Press, New York 
(1969). 

S. Brazovskii, J.-P. Pouget, S. Ravy and S. Rouziere, 
Phys. Rev. B55, 3426 (1997). 

S. Rouziere, S. Ravy, S. Brazovskii, and J.-P. Pouget, in 
preparation. 

For a review see a) C. Schlenker, Low-Dimensional Elec- 
tronic Properties of Molybdenum Bronzes and Oxides 
(Kluwer Academic Publishers, Dordrecht, 1989), p. 159 
and b) J.-P. Pouget, ibidem, p. 87. 

J.-P Pouget, C. Noguera, A.H. Moudden and R. Moret, 
J. Phys. (France) 46, 1731 (1985). 

S. Rouziere, S. Ravy and J.-P. Pouget, Synth. Metals 86, 
2131 (1997). 

Due to the lack of a model for fitting the whole profile, 
we chose the simplest way to evaluate £&* . In Rouziere 
et al, Synth. Metals 70, 1259 (1995), the Lorentzian- 
squared profile was used (£ = 0.6436/ HW H M). 
J. M. Ziman, Principles of the Theory of Solids, Cam- 
bridge University Press 1972. 

I. Tutto and A. Zawadovskii, Phys. Rev. B32, 2449 
(1985); J. R. Tucker, Phys. Rev B40, 5447 (1989). 



4 



